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In granular superconductors, individual grains can contain bound Cooper pairs while the system 
as a whole is strongly insulating. In such cases the conductivity is determined by electron hopping 
between localized states in individual grains. Here we examine a model of hopping conductivity in 
such an insulating granular superconductor, where disorder is assumed to be provided by random 
charges embedded in the insulating gaps between grains. We use computer simulations to calculate 
the single-electron and electron pair density of states at different values of the superconducting 
gap A, and we identify "triptych" symmetries and scaling relations between them. At a particular 
critical value of A, one can define an effective charge V2e that characterizes the density of states 
and the hopping transport. We discuss the implications of our results for magnetoresistance and 
tunneling experiments. 



I. INTRODUCTION 

Granular superconductors are arrays of superconduct- 
ing granules that are connected by electron tunneling. As 
such, these systems combine the unique electronic spec- 
trum of superconducting quantum dots with the strong 
Coulomb correlations that are ubiquitous in disordered 
systems 1 . Among the more celebrated properties of 
granular superconductors are a giant magnetoresistance 
pealP^ and a superconductor-insulator transition that 
can be tuned by disorder or magnetic fields. So far, a 
comprehensive theory of the electron conductivity that 
can explain these features remains elusive. 

In this paper, we focus on the strongly disordered 
limit, where the array of superconducting grains as a 
whole is insulating while individual grains may still retain 
prominent features of superconductivity. In this case, 
electronic states are localized and electron conduction 
proceeds by phonon-assisted tunneling, or "hopping," 
of electrons between grains through the insulating gaps 
which separate them. In principle, electronic conduction 
can occur either through tunneling of single electrons or 
through simultaneous tunneling of an electron pair. In 
this introductory discussion we concentrate primarily on 
hopping by single electrons; hopping by electron pairs is 
discussed more fully in the body of the text. Here we note 
only that coherent tunneling of Cooper pairs (the Joseph- 
son effect) is neglected throughout this paper, since it is 
not relevant in the strongly disordered limit that we are 
considering. 

Since hopping conductivity is a thermally activated 
process, its magnitude at a given temperature T depends 
on two important energy scales associated with the spec- 
trum of electron energy states within each grain. The 
first is the charging energy E c = e 2 /2Co, where e is the 
electron charge and Co is the self-capacitance of a sin- 
gle grain. The importance of the charging energy can be 
seen by considering that, in a neutral system, conduction 
requires an electron to hop from one neutral grain to an- 
other, thereby producing two charged grains, each with 



Coulomb self-energy E c . The second important energy 
scale is the superconducting gap A, which represents an 
activation energy for separating a Cooper pair. In the 
limit whe re A /E c 0, the array is equivalent to a gran- 
ular metaPEHU i n the opposite limit, where A/E c oo, 
each grain has the properties of a bulk superconductor. 
In this paper our focus is on exploring the novel physics 
that results when E c and A are similar in magnitude. 

Since the superconduct ing gap A is typically on the 
order of 1 meV or smaller^^, E c ~ A implies that the 
self-capacitance Co > 80 aF. This relatively large self- 
capacitance can be achieved either by fabricating large 
grains or by surrounding the grains by an environment 
with a high effective dielectric constant ft, so that the 
product of ft and the grain diameter D satisfies kD > 400 
nm. For three-dimensional (3d) arrays, large Co can also 
be achieved by making an array of very densely-packed 
grains, for example, cubic grains separated by a thin in- 
sulating layeJ^. In this paper we assume that the Joseph- 
son coupling energy J between grains satisfies J <C E c , 
so that the array is indeed insulating 1 -! regardless of the 
value of A, and coherent tunneling of Cooper pairs is ab- 
sent. Since we are considering the case of relatively large 
grains, we also assume that the spacing S between dis- 
crete electron energy eigenstates within the grain satisfies 
6 < E c . 

In relatively dense arrays, electron conduction can oc- 
cur both through hopping of electrons between nearest- 
neighboring grains and through "cotunneling" of elec- 
trons between distant grains via a chain of intermediate 
virtual states in intervening grainj ^ n l 12 l In the pres- 
ence of some disorder, the latter mechanism dominates at 
low temperatures, where the length of hops grows to opti- 
mize the conductivity. This transport mechanism, called 
variable range hopping (VRH), was introduced by MotlP 1 
and is ubiquitous in low-temperature systems with local- 
ized electron states. 

The optimum average hop length at a given tempera- 
ture depends on the distribution of available ground state 
energies among grains in the system, called the "density 
of ground states" (DOGS). For systems with unscreened 
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Coulomb interactions, it can be shown on the basis of a 
very general stability criterion of the global ground state 
that the DOGS must vanish^ at the Fermi level [i. This 
vanishing DOGS is called the Coulomb gap, and in the 
canonical Coulomb glass model of disordered systems it 
leads to a conductivity a that obeys the Efros-Shklovskii 
(ES) law: 

acxexpHTEs/T) 1 ^ (1) 

where 



is a characteristic temperature, C is a numerical coef- 
ficient, and £ is the electron localization length of the 
array. Eq. ([!]) has been observed in a number of granular 
metals and superconductors at low temperature (see Ref. 
[Hand references therein). 

In a previous papei!^, we used a computer simulation 
to explore VRH in two-dimensional (2d) and 3d arrays 
of monodisperse normal metallic grains. In such systems 
disorder is provided by donors and acceptors randomly 
situated in the interstitial spaces between grains — for ex- 
ample, in the metal oxide of the grains. We showed that 
as a consequence of the periodic charging spectrum of 
individual grains there is not one but three identical ad- 
jacent Coulomb gaps in the DOGS (one full gap at the 
Fermi level and two "half-gaps" on either side) , which to- 
gether form a structure that we termed a "Coulomb gap 
triptych." Unlike in conventional Coulomb glass models, 
in metallic granular arrays the DOGS has a fixed width 
in the limit of large disorder. 

This previous study can be considered as a model for 
a granular superconductor in the limit where A/E c — > 0. 
In the present paper, we generalize the theory of Ref. 9 
to the case of finite A* = A/E c . Specifically, we assume 
that within each grain electrons can form Cooper pairs, 
thereby lowering the system energy by — 2A per pair. We 
use this model to study the DOGS and conductivity as a 
function of the gap A and the temperature T . Our result 
for the DOGS of single electrons, which we denote gi(E), 
is shown in Fig. [T^i for different values of A*. Fig. [TJd 
shows the DOGS for electron pairs, gi{E). In both cases 
the dimensionless energy e = E/E c is defined relative to 
the Fermi level, so that states with e < are occupied 
and states with e > are empty. These results for the 
DOGS are explained more completely below, but here we 
note briefly that for sufficiently large disorder the system 
exhibits repeated Coulomb gaps at all values of A, and 
the width of the DOGS is determined by the values of 
E c and A and not by the strength of the disorder. 

For energies close to the Fermi level our results for 
gi(E) and g2(E) are similar to those of an earlier semi- 
nal worlP^, which aimed to capture the effect of pairwise 
attraction of electrons on the Coulomb gap and VRH 
conductivity. The authors of Ref. [TS] started from the 
canonical Efros model of the Coulomb glass^with strong 
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FIG. 1. (Color online) Single electron and pair DOGS, g*(e) 
and 32(e), of a regular 2d array of monodisperse metallic 
grains as a function of the dimensionless electron energy 
e — E / E c at different values of the superconducting gap 
A* = A/E c . At A* < 1, the single electron DOGS g* has 
a soft Coulomb gap at e — 0, while the pair DOGS has a 
hard gap, and the situation is reversed for A* > 1. A* = I is 
a critical point at which both g\ 2 have a soft Coulomb gap. 
The three DOGS curves corresponding to g*(e) at A* = 0, 1 
and 32(e) at A* > I constitute "Coulomb gap triptychs" and 
can be scaled onto each other by rescaling the electron charge, 
as discussed in Sec. |III| One can equivalently say that these 
three curves exhibit effective charges e, \/2e, and 2e, respec- 
tively. 

disorder and added the possibility of occupation of a site 
by two electrons with a finite (positive or negative) in- 
teraction energy U. They used this model to study how 
varying the on-site energy U affects gi(E), g2(E), and 
the hopping conductivity a in the presence of large ex- 
ternal disorder. In the present paper we examine a model 
that is more realistic for granular superconductors, and 
we confirm a number of interesting observations made in 
Ref. [H 

The most interesting of these observations is the ap- 
pearance of effective charges e* = fe, e* = 2e, and 
e* = V2e at A* < I, A* > 1, and A* = 1, respec- 
tively. These effective charges play a prominent role in 
the DOGS and the hopping conductivity, as we show be- 
low. The unusual V2e charge that arises at A* = I is 
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the product of a degeneracy in the electronic spectrum, 
and can be seen as a consequence of singly-charged elec- 
trons hopping in a Coulomb landscape whose properties 
are determined by doubly-charged Cooper pairs^. 

The remainder of this paper is organized as follows. 
In Sec. [Il] we define the system being studied and out- 
line our simulation technique. In Sec. |III| we describe our 
main results for the DOGS and conductivity, focusing 
primarily on 2d arrays, and we present their implications 
for magnetoresistance. In Sec. |IV| we show that our re- 
sults generalize to the 3d case as well. Sec. [V] is devoted 
to translating our results for the DOGS into a predic- 
tion for tunneling experiments. We close in Sec. VI with 
concluding remarks. 



II. MODEL 

In this paper we consider an array of identical, spher- 
ical grains with diameter D arranged in a regular, d- 
dimensional square lattice with lattice constant D' > D. 
For simplicity of discussion, during the majority of this 
paper we focus on case d — 2; results for d = 3 are pre- 



sented in Sec. IV Disorder in this system is assumed to 
be provided by impurity charges ±e that are embedded 
in the insulating interstitial spaces between grains. Such 
impurity charges can be thought to effectively create a 
random fractional donor charge Qi that resides on each 
grain i, for reasons that are explained more fully in Ref. 
9; The net charge of the grain can then be written as 
Qi = Qi~ eni, where n, is the integer number of electrons 
that reside on the grain relative to its neutral state. We 
emphasize that rii can be a positive or negative integer, 
and can be defined as n, = Ni — Ii, where Ii is the num- 
ber of positive ions and Ni the number of electrons at 
grain i. Within each grain, the Ni electrons can form 
bound pairs through the local attraction energy A. In 
general, A can be tuned by an applied magnetic field B, 
as discussed below. 

Given this model, the Hamiltonian for the system can 
be written 



(Qi - eriif 



2C 



^CT-.^Qi-em^Qj-erij) (3) 

(id) 



Ni 
2 



Here, the first term describes the Coulomb self-energy of 
each grain and the second term describes the Coulomb 
interaction between charged grains. The coefficient Cy 
is the inverse of the matrix of electrostatic induction Cy . 
In the numerical simulations that we describe below, we 
make the approximations Cq — kD/2 and C^ 1 = l/fcry. 
The third term in the Hamiltonian describes the total 
pairing energy of electrons; Ni is the number of electrons 
and |_7V^ / 2J is the number of electron pairs within grain 
i. In Ref. 1151 the authors proposed a similar Hamiltonian 
as a model for disordered superconducting films such as 



InOa;. Unlike in Eq. Q, the model considered in Ref. 
1151 assumes that the electron occupation numbers are re- 
stricted to Ni = 0, 1, 2 and that disorder is provided by 
random, uncorrelated site energies rather than by the 
random charges Qi. While we consider our model more 
realistic for granular superconductors, we will show that 
it reproduces many of the features reported in Ref. 1151 

Because of the presence of the impurity charges, elec- 
trons become redistributed among grains from their neu- 
tral state in order to screen the disorder Coulomb po- 
tential. The corresponding ground state arrangement of 
electrons among grains plays an essential role in the con- 
ductivity, since it determines the lowest empty and high- 
est filled electron energy levels at each grain. In our 
numerical simulation described below, we search for the 
set {rii} that minimizes the Hamiltonian and use it to 
calculate the DOGS and the conductivity. 

In conventional Coulomb glass models, the character- 
istic strength of the disorder is a free p arameter that 
determines the width of the DOGS^ES. One can ex- 
pect that in our problem a similar role is played by the 
typical magnitude of the disorder charge Qi, which re- 
flects the average number of impurity charges surround- 
ing each grain. In fact, however, in the limit where there 
are many such charges one can effectively adopt a sim- 
ple model in which the value of Qi is chosen randomly 
from the uniform distribution Qi 6 [— e, +e\. To see why 
this model is valid, consider that each grain minimizes 
its Coulomb self-energy by minimizing the magnitude of 
its net charge, |<jfj| = |Qj — enj|. In the absence of any 
Cooper pairing, m may freely take any integer value in 
order to arrive at a state for which — e/2 < qi < e/2. 
If one assumes, on the other hand, that Cooper pairing 
is so strong that all electrons are paired in the ground 
state (A* 3> 1) , then Ni = rii + Ii may still freely 
take any even-integered value, so that in the ground state 
— e < qi < e. In other words, regardless of the value of A, 
each grain can effectively adjust to the presence of an ar- 
bitrarily strong charge disorder by changing its electron 
number n, so that its net charge acquires a magnitude 
smaller than e. For this reason, in the limit of large disor- 
der the DOGS has a fixed width, as first explained in Ref. 
HI For the results presented below, we take Qi to be ran- 
domly chosen from the uniform distribution Qi G [— e, e]. 
The ion number Ii is assumed to be very large, so that 
electrons are never completely depleted from any given 
grain. Zj is also taken to be even or odd with equal prob- 
ability; the relevance of this choice is explained below. 

In our analysis below it is convenient to introduce the 
following dimensionless units, which reduce the number 
of free variables in the problem. We introduce the di- 
mensionless distance between the centers of grains i and 
h 



D ' 



the dimensionless charge 



Q l 



(4) 



(5) 
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the dimensionless electron energy 

e = E/E c , 



(6) 



the dimensionless DOGS for single electrons and electron 
pairs 



g*{e) = E c D d g h2 (e), 



the dimensionless temperature 



T* 



2Dk B T 



and the dimensionless resistivity 



(7) 



(8) 



(9) 



where po is a prefactor for the resistivity with a weak, 
power-law dependence on temperature. We also assume 
that the gap between neighboring grains D' — D -C D, 
so that D' ~ D. The problem then loses any explicit 
dependence on the diameter or the localization length. 
With these definitions, one can write the Hamiltonian of 
Eq. Q in dimensionless units as 



(id) ij 



E 
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(10) 



If one is given the ground state electron occupation 
numbers {rii}, then one can determine the highest oc- 
cupied electron energy state, e\~ , and the lowest empty 
state, e\ + , at a given grain i. These energies determine 
the contribution of the grain i to the single-electron con- 
ductivity, and are given by: 



i 





2A* 



2A* 




A, odd 
JVi even 



iVj odd 
iVi even 



(11) 



(12) 



From Eqs. (Ill and (12) one can see that the spec- 
trum of single-electron energy levels at a given grain i 
depends on the "parity" of the grain: whether the num- 
ber of electrons in the neutral state, Jj, is odd or even. 
To understand why this is the case, consider first the 
spectrum of a single grain with A = 0. In such a grain, 
e\ + — e\~ = 2 regardless of the number of electrons in 
the grain. This implies a ladder of electron energy levels, 
spaced by 2E C , corresponding to different charge states 
of the grain. These energy levels are shown schematically 
in the left side of Fig. [2j When A* is finite, on the other 
hand, those energy states corresponding to an even to- 
tal number of electrons in the grain become shifted by 
— 2A* as a consequence of the attractive interaction be- 
tween electron pairs. As a result, e 1+ — e\~ = 2 — 2A* 



for grains with odd A, and e\ + — e\~~ = 2 + 2A* for 
grains with even A,. This suggests that the energy to 
add or remove one electron from the grain's neutral state 
depends on the parity of the grain, as shown in the cen- 
ter of Fig. [2] (The importance of the grain parity for 
its electronic sp ectrum has been well e stabli shed by pre- 
vious theoretica l 18 * 19 ^ and experimenta l 20 * 21 ! studies.) At 
A* = 1, pairs of electron energy states become two-fold 
degenerate, as shown on the right side of Fig. [2] As a 
consequence, at A* > 1 in the ground state all grains 
have an even total number of electrons, regardless of the 
disorder strength. This uniform pairing has an important 
consequence for the DOGS, as discussed below. 
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FIG. 2. Single-electron energy levels of an isolated neutral 
grain. At A* = (left), the Coulomb self-energy produces 
a spectrum where different charge states e are separated in 
energy by 2. At < A* < 1 (center), those energy levels 
corresponding to the addition of an electron to a grain with 
an odd total number of electrons are shifted by — 2A*. At 
A* = 1 (right), the difference in self-energy 2E C between two 
successive charge states is compensated by the pairing energy 
— 2A, so that pairs of subsequent electron levels merge. The 
Fermi level at e — is indicated schematically by a dashed 
line. 

The diagram of Fig. [2] shows that the spacing between 
energy levels at A* = 1 is doubled relative to that of 
A* = 0. One can observe that this same increased spac- 
ing could be achieved if the electron charge e were re- 
placed with an effective charge V2e, so that the charging 
energy E c oc e 2 is doubled. In fact, this effective charge 
y/2e plays a prominent role for the DOGS and conduc- 
tivity at A* = 1, as will be shown in Sec. |III| 

In the presence of some disorder, the ladder of energy 
states depicted in Fig. [2] becomes shifted randomly up or 
down from one grain to the next by the disorder poten- 
tial. The values of e\~ and e\ + for each grain — those en- 
ergy states just above and just below the constant global 
Fermi level — contribute to the DOGS g*{e). 

Thus far we have focused our discussion on hopping 
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by single electrons, which is characterized by a localiza- 
tion length £ = £1. In principle, conduction may occur 
through simultaneous hopping of an electron pair as well, 
with a distinct localization length £ = £2- For example, 
one can expect pair tunneling to become dominant in 
the limit A* — > 00, where thermally activated breaking 
of bound Cooper pairs is completely suppressed. In or- 
der to discuss conduction by electron pairs, one can sim- 
ilarly define the energy associated with pair excitations, 
in analogy with Eqs. (12) and (11). Specifically: 



= -4tf-4-2££-2A 



.2+ 



= -4 9 * + 4-2^-%-2A*. 



r ' 



(13) 
(14) 



Note that, unlike for single electron excitations, for pairs 
we have £ 2+ — e 2- — 8 regardless of the parity of the 
grain. This suggests that the ladder of energy states cor- 
responding to pair excitations has a uniform spacing 8E C , 
and thus all pair excitation energies are independent of 
A. This is as expected, since the total number of bound 
pairs in the system is unchanged by the simultaneous 
tunneling of a pair. As with the single electron energy 
levels, the disorder potential produces a random shifting 
of the two-electron energy levels from one grain to an- 
other. The energies e 2- and ef + in the ground state are 
histogrammed to produce the pair DOGS, g% (e). 

In order to evaluate numerically the DOGS, we use a 
computer simulation to search for the ground state ar- 
rangement of electrons, {n^}, in a finite array of grains. 
This is done by looping over all pairs ij and attempt- 
ing to move either one or two electrons from i to j. If 
the proposed move lowers the total system energy H* , 
then it is accepted, otherwise it is rejected. This process 
is continued until no single-electron or pair transfers are 
possible that lower H*. Equivalently, one can say that 
for all i, j we check that two sets of ES ground state 
criteria are satisfied: 



and 



— £■ 



— E s 



1/r*, > 0. 



- 4/r* ? > 0. 



(15) 



(16) 



The final arrangement of electrons can be called a 
"pseudo-ground state," which is not strictly equal to the 
true ground state of the system but which general ly pro- 
vides an identical DOGS up to very small energies^- 17 * 22 * 23 !. 

Once the energies {e\ ' 2 } are known, we evaluate the 
resistivity using the approach of the Miller-Abrahams 
resistor networlPH. This approach is described in de- 
tail in Refs. OH and [351 but here we give a brief con- 
ceptual overview. In the Miller-Abrahams description, 
each pair of grains ij is said to be connected by some 
equivalent resistance R^. The value of Rij increases ex- 
ponentially with the distance between the grains and 



with the activation energy Eij required for hopping be- 
tween them according to Rij cx exp[2r,j/£ + Eij/ksT]. 
Note that, using the dimensionless units of Eqs. Q - 
([9]), one can define the dimensionless logarithm of the 
resistance lni?*,- = r*j +Eij/T* , which has no explicit de- 
pendence on £. The resistivity of the system as a whole is 
found using a percolation approach. Specifically, we find 
the minimum value R c such that if all resistances Rij 
with Rij < R c are left intact while others are eliminated 
(replaced with R = 00), then there exists a percolation 
pathway connecting opposite faces of the simulation vol- 
ume. The system resistivity is equated with R c D d ~ 2 . 

In principle, single-electron hopping and pair hopping 
provide parallel mechanisms for charge transport be- 
tween a given pair of grains ij, and so they can be rep- 
resented as parallel resistors connecting the two grains. 
In most situations, however, one of the two mechanisms 
dominates the conductivity while the other can be ne- 
glected, as we show below. We therefore focus primar- 
ily on the case where single and pair excitations can be 
treated as independent, non-connected resistor networks 
with resistivities p\ and P2, respectively. Some limited 
results for mixed conduction are provided at the end of 
the following section. 

All numerical results for 2d systems presented in the 
following section correspond to simulations of 100 x 100 
lattice sites with open boundary conditions, averaged 
over 1000 independent, random realizations of the dis- 
order. Energies are defined relative to the Fermi level 
(j,, so that in the ground state e\~ < 0, < and 
e\ + > 0, e 2+ > for all i. 



III. RESULTS 

The DOGS is shown in Fig. [I] for single electron exci- 
tations, gl(s), and for pair excitations, 32( e )i a ^ different 
values of the gap A*. For each curve, the DOGS vanishes 
at the Fermi level (e — 0), as required by the stability 
criteria of Eqs. ( 15 ) and ( 16 ). One can also note that the 
DOGS generally becomes wider with increasing A* as a 



consequence of the widening gaps between odd and even 
electron energy levels (see Fig. [2| . The evolution of the 
DOGS with A* can be understood more completely as 
follows. 

At A* = 0, the array is equivalent to a granular normal 
metal. As a consequence, the curve g\{e) at A* = is 
identical to the one reported in Ref. El The most salient 
feature of this curve is its "triptych" symmetry, with two 
identical peaks that are symmetric about their centers. 
As explained in Ref. '9", this symmetry is a result of the ES 



stability criterion of Eq. ( 15 ) in conjunction with the uni- 
form spacing between electron energy levels (at A* = 0, 
e\ + — e\~ = 2 for all i, as shown in Fig. [2]). Thus, the 
"soft" Coulomb gap at e = gets repeated identically 
at e = ±2. On the other hand, the pair DOGS, g^i?), 
at A* = has a hard gap at the Fermi level with width 
4. This hard gap can be understood by considering that 



() 



at A* = 0, Eqs. (11) - (14| imply that e\ 



_2± 



24± ± 2. 

Since e\ < and e\ + > for all i, we have |e 2± | > 2, 
and therefore there must be a hard gap of width 4. Phys- 
ically, one can say that the gap arises in g\ (e) because the 
charging energy 4E C associated with adding two electrons 
to a given grain is larger in magnitude than the random 
Coulomb potential, which is screened effectively by the 
rearrangement of single electrons. As a consequence of 
the relation between e? and e\ , at A* = the two 
DOGS can be mapped onto each other via the relation 
gi(e) = 2<72 [2e + 2sgn(e)]. A slightly different version of 
this relation was reported in Ref. [15] 

When the pairing interaction is finite but small, < 
A* < 1, gl(e) is unchanged very close to the Fermi 
level, but away from the Fermi level it becomes some- 
what broadened due to the widening energy gaps between 
even and odd parity electron states (see Fig.[2]). The pair 
DOGS, meanwhile, retains a hard gap near the Fermi 
level, but the width of this gap shrinks to 4(1 — A*). 

In the opposite case, where the pairing interaction is 
strong enough that A* > 1, the situation is reversed. 
That is, the single-electron DOGS g*(s) acquires a hard 
gap at the Fermi level while g% (e) has only a soft Coulomb 
gap. This result can be understood by first noting that 
at A* > 1, all electrons are paired in the ground state. 
This is true because at A* > 1 any grain with an odd 
number of electrons can lower its energy by acquiring an 
electron from a distant grain with electron energy close 
to the Fermi level (or from the voltage source). Making 
use of Eqs. (11) - (14) for even-parity grains produces 



the relation £„■ = 



if* ± (A* - 1). Since e-+ > and 



e?~ < 0, it follows that all Ie]*] > A 



1 for all i, and 
thus there is a hard gap in g\{e) of width 2(A* — 1). 
Physically, this hard gap arises because the pairing in- 
teraction is stronger than the disorder Coulomb poten- 
tial, which is screened effectively by Cooper pairs. Thus, 
any excitation of single-electron hops requires a finite ac- 
tivation energy of at least A* — 1. The relations be- 
tween £ 2± and at A* > 1 imply a mapping between 
g*(s) and <?2 (e) that was also noticed by Ref. [T51 namely 
52( e ) = |ff*[| £ + sgn(e)(A* - 1)]. At such large values 
of A*, the fixed relation e 2+ — ef~ = 8 implies that g\ (e) 
becomes saturated and has a fixed width for all A* > 1. 

At the point where A* = 1 precisely, some remarkable 
features emerge in the DOGS. This might be expected by 
noticing the special role played by A* = 1 in the single- 
electron energy spectrum; this is the point where pairs of 
energy levels become degenerate (see Fig. [2j right). At 
A* = 1 neither g\{e) nor g% (e) has a hard gap, and in 
fact the two DOGS can be mapped onto each other via 
the simple relation g\ (e) = \g\{\e). In addition, there 
is a simple scaling relation between g\{e) at A* = 1 and 
g\{e) at A* = 0. Namely, 



1 



;9t (e/2) 



= 2ff2(2£) 



A*>1 



(17) 



The second equality in Eq. (17) can be understood 
in a straightforward way. Indeed, the second equality 



suggests that one can arrive at the pair DOGS at large 
A* by taking the single-electron DOGS at A* = and 
rescaling the value of the electron charge by a factor of 
2. Replacing e by an effective charge e* = 2e in the 
unit of energy E c produces a factor 4 expansion of the 
x-axis and a factor 4 contraction of the y-axis, which 
is equivalent to the second equality in Eq. (17). This 



scaling can be expected, since for large pairing interaction 
A* > 1, all electrons are paired, and one can naturally 
think that only charge 2e objects exist in the problem. 
Thus, at such large A* the problem of the arrangement 
of electron pairs in the disorder potential is equivalent to 
the problem of the arrangement of single electrons in a 
disorder potential, with rescaled units. 



The first equality in Eq. (17), on the other hand, is un- 
expected, since it implies that g\{e) | A ,_ 1 can be deter- 
mined from gl(s) | A „_ by replacing the electron charge 

with an effective charge e* = y/2e. This remarkable fea- 
ture of the DOGS at A* = 1 was first pointed out by 
Ref. HH1 Those authors showed that the result e* = \/2e 
is the natural consequence of single electrons hopping 
in a Coulomb landscape that is shaped predominantly 
by Cooper pairs. More formally, one can say that the 
pair stability criterion of Eq. ( fl6| ) produces a stronger 
cons trai nt on g\{e) than the single-electron criterion of 



Eq. qisp 

whic 



This can be seen by substituting e? = 2e 



_i± 



1 is correct at A* = 1 (see Fig. v2b 



1+ 



As a result, one finds that e 
dimensionfull units, £7- — E) 



2/r- 



into Eq. (16) 
> 0, or ii 



-3 -1 « 
^ - {V2e) 2 /K riJ > 0. Re- 
peating the traditional derivation of the Coulomb gap ^ 
starting with this inequality leads to an effective charge 
e* = \/2e in the DOGS. 

In addition to its importance for the DOGS, the effec- 
tive charge e* plays a prominent role in the electron con- 
ductivity. Specifically, it enters the characteristic tem- 
perature T E s in the ES law [Eq. Q]. Since T E s oc e 2 [see 
Eq. pb], the arguments above suggest that if one defines 
the ES temperature Tg S (A*) for single-electron conduc- 
tivity and Tg S (A*) for pair conductivity at a given value 
of A* , then these should satisfy 

Tg s (A* > 1) = 2T| S (A* = 1) = 4T! S (A* = 0). (18) 

In order to verify this prediction, we measured the 
single-electron resistivity In p* and the electron pair re- 
sistivity In P2 at various values of A* and over a range 
of temperatures using the resistor network approach de- 
scribed in Sec. [TT1 The result is plotted in Fig. [3] as 
\np* versus (T*) / 2 . As expected, at low temperature, 
T* <C 1, the conductivity is well described by the ES 
law in all cases. By making linear best fits to the data 
at low temperature, we find that the corresponding tem- 
peratures T E s indeed satisfy Eq. (18). This can be seen 



from the dashed lines in Fig. [3j which show three fit lines 
with relative slopes 1 : y2 : 2, as predicted by the cor- 
responding effective charges e*. If the data in Fig. [3] 
are fitted with independent best fit lines, we find that 



T| S (A* > 1) 



2.2T! S (A* 



4.3T| S (A* 



0), 
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which is within our numerical uncertainty of the predic- 
tion in Eq. (fl8b . 



Tes = ^es(O)- According to the second equality in Eq. 
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FIG. 3. (Color online) The temperature dependence of the 
resistivity for single electron conduction at A* = 0, 1 and pair 
conduction at A* > 1. The dimensionless resistance hip* 
is plotted against (T*) -1 / 2 to illustrate that the resistivity 
follows the ES law [Eq. Q] at low temperatures. The dashed 
lines are linear fits whose slopes have the ratio l:y/2:2. 

The evolution of Tes with A suggests an interest- 
ing mechanism for the magnetoresistance of the sample. 
Generally speaking, the pairing energy A in a supercon- 
ducting material decreases monotonically 2 ^ with the in- 
tensity of an applied magnetic field B. Thus, by ap- 
plying a magnetic field one can tune the pairing energy 
and thereby alter the DOGS, the ES temperature Tes, 
and the resistivity. In the following discussion we assume 
that this tuning of A is the primary role of an applied 
magnetic field, and we ignore the effect of the magnetic 
field on hopping interference phenomena^. One could 
also imagine that the magnetic field is applied parallel to 
the array, so that all hopping trajectories encircle zero 
magnetic flux. 

In order to investigate this mechanism for magnetore- 
sistance, we consider first the case where all conduction 
is due to single electron hopping. This would be the 
case, for example, when £,2/^1 "C 1. In such a case the 
results of Figs, [l] and [3] imply a monotonic negative mag- 
netoresistance. That is, as a magnetic field B is applied, 
the gap A decreases, leading to a larger DOGS near the 
Fermi level and thus to enhanced conductivity. More 
specifically, if the superconducting gap is large enough 
that at zero magnetic field A* > 1, then in the ab- 
sence of a magnetic field the single-electron DOGS has 
a hard gap. This hard gap implies that at low temper- 
atures T* <C (A* — 1), the resistivity is very large and 
described by an Arrhenius-type activation law. When 
B is increased to the point that A* = 1, the resistiv- 
ity becomes smaller and obeys the ES law with a char- 
acteristic temperature Tg s (l). As the magnetic field 
is increased even further, Tes decreases and the resis- 
tivity declines. This decline continues until such large 
fields that A* « 1, when the resistivity plateaus and 



(18), at small temperatures one should expect that the 
large- B resistivity and the resistivity at A* = 1 are re- 
lated by [lnp*(A* = l)]/[ln^(A* = 0)] ~ y/2. 

This result is confirmed in Fig. [4j which shows the 
single-electron resistivity as a function of the supercon- 
ducting gap A* at various values of temperature. As 
expected, the resistivity indeed declines with decreasing 
gap (increasing B), and at very small temperatures (large 
(T*)- 1 / 2 ) the relation [lnpJ(A* = l)]/[lnp*(A* = 0)] = 
y/2 is nearly satisfied. This result provides an additional 
confirmation of the picture of an effective electron charge 
y/2e at A* = 1. 
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FIG. 4. (Color online) Resistivity for single-electron hopping, 
In pi , as a function of the superconducting gap A* at different 
values of the temperature T* . The resistivity is normalized 
by its value at A* = 0, and one can see that for small tem- 
peratures the ratio [In pi (A* = l)]/[lnp*(A* = 0)] seems to 
approach y/2. The declining resistivity with decreasing gap 
implies a negative magnetoresistance. The dotted vertical line 
indicates A* = 1, which can be thought of as the point where 
the resistivity crosses over from an activated dependence to 
the ES law with increasing magnetic field (at small tempera- 
ture). 

We would like to emphasize that this mechanism for 
negative magnetoresistance is quite unusual, and cannot 
be understood simply as a reduction of some activation 
energy due to weaker Cooper pairing. Rather, the nega- 
tive magnetoresistance arises because decreased A* leads 
to a DOGS g\{e) that is less depleted by intimidation by 
Cooper pairs, and thus to enhanced electron conduction 
at low temperature. 

The results of Fig. [4] focus on the case where conduc- 
tion is provided by single electrons only, which is ap- 
propriate when £2/^1 *C 1. On the other hand, when 
the localization lengths £1 and £2 a re similar in mag- 
nitude, the conduction should be dominated by single- 
electron hopping at A* -C 1 and by pair hopping at 
A* 1. This is the case because at all A* ^ 1 one 
of the two DOGS has a hard gap. By increasing the 
magnetic field, then, one can apparently produce a tran- 



sition between pair-dominated conduction and single- 
electron-dominated conduction, provided that A* > 1 
in the absence of applied field. Such a transition may 
help to explain the giant magnetoresistance peak seen in 
experiments^"^, as was proposed by Ref. 1151 

To investigate this possibility, we performed simula- 
tions to measure the resistivity at different values of the 
localization lengths £1, £2 and the temperature T, using 
a resistor network that allows for mixed conductivity of 
singles and pairs. We find that if £1, £2> and T are chosen 
such that the resistivities are nearly equal at A* ^S> 1 and 
A* = 0, then one can indeed observe a moderate peak in 
the resistivity in the vicinity of A* = 1. One such result 
is shown in Fig. [5j and is qualitatively similar to a result 
obtained in Ref. JTS] for an array in which the quantity 
E c — A varies strongly between grains. 
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FIG. 5. (Color online) The resistivity, p/po, as a function of 
the superconducting gap A* for a 2d array with localization 
lengths £1 = D and £2 = WD and temperature T = 0.1E c /1cb 
[so that (Ti)~ 1/2 = 2.1 and (T 2 * )" 1/2 = 6.7]. The maximum 
in p/po suggests a magnetoresistance peak associated with 
the transition from pair-dominated conduction (at large A*, 
small magnetic field) to single electron-dominated conduction 
(at small A*, large magnetic field). 



While this result is promising, we caution that by itself 
it does not provide a satisfactory qualitative description 
of the magnetoresistance peak observed in experiment. 
For example, the peak in Fig. [5] arises out of the deeply 
insulating state, p ~ I0 3 po. Since the constant po is 
generally on the order of h/e 2 w 26 kfi, this disagrees 
with experimentPSl, where the magnetoresistance peak is 
seen to arise from a state with p ~ h/e 2 . The appearance 
of a noticeable peak also apparently requires a large ratio 
which is likely to be possible only very close to the 
superconductor-insulator transition. Such large values of 
£2 probably go beyond the limit of applicability of our 
model. 



IV. 3D ARRAYS 

Thus far our presentation of results has focused on the 
case of 2d arrays. In this section we briefly report on 
simulations of the DOGS and resistivity in 3d arrays. 
Generally speaking, while some details of the shape of 
the DOGS and the magnitude of the resistivity are mod- 
ified relative to the 2d case, the triptych structure of 
the DOGS and the values of the effective charges remain 
unchanged. All results in this section correspond to sim- 
ulated 3d systems of 24 x 24 x 24 lattice sites with open 
boundary conditions, averaged over 1000 realizations of 
the disorder. 

When considering the DOGS, the most prominent dif- 
ference between 2d and 3d systems is that in 3d the 
ES criterion [Eq. ( 15 )] imposes a stronger constraint 



on g\ ( e). Specifically, in d dimensions the ES criterion 
implies^ that g*(e) < fcd|e| d_1 , where kd is a constant, 
so that in 3d the DOGS vanishes at least quadratically 
with energy near the Fermi level while in 2d it is con- 
strained to vanish only linearly. Fig. [6] shows gl (e) and 
52(e) for 3d arrays, and one can see that in cases where 
<?J (e) is ungapped it indeed vanishes as a higher power of 
e near the Fermi level. 

Nonetheless, the most important qualitative features of 
the DOGS from the 2d case remain in 3d as well. Specifi- 
cally, the curves corresponding to g\(e) at A* = 0, 1 and 
52(e) at A* > 1 have the "triptych" structure of two 
identical, symmetric peaks, and they can be scaled onto 



each other using the same scaling relations of Eq. (17 1 



This implies that in 3d we have the same effective charges 
e* = le, e* = \/2e, and e* = 2e for single electrons at 
A* = 0, 1 and for pairs at A* > 1, respectively. 

Thus, the most important conclusion from our results 
in 2d remains for the 3d case as well. This is as expected, 
since, as explained in Sec. |III| the effective charges arise 
from the single-electron energy spectrum (Fig. [2]), and 
are therefore independent of the system dimensionality. 

One can also check that in 3d the effective charges 
have the same influence on the ES temperature as pre- 
dicted by Eq. (18 1. By numerically evaluating the resis- 
tivity of these 3d systems, we indeed find that the ES 
temperatures Tgs obey Eq. (fl8|). A plot of the dimen- 



sionless resistivity \np* against (T*) -1 / 2 in 3d is essen- 
tially identical to that of Fig. [3j with slight downward 
shifts in the magnitude of the resistivity relative to the 
2d case. Making independent linear fits to the data gives 
T|g(A* > 1) « 2.6T!JA* = 1) « 5.1T! S (A* = 0), 
which agrees with Eq. ( |18[ ) to within our numerical un- 
certainty. 



V. TUNNELING EXPERIMENTS 

In the previous sections we presented results for the 
DOGS and we showed that these results have important 
consequences for the characteristic temperature Tes and 
for the magnetoresistance. In this section we discuss how 
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Since the spacing d between discrete electron energy 
levels within the grain satisfies 6 -C E c , as explained 
in the Introduction, we can take the density of states 
f(E) within each grain to be a continuous function. For 
metallic grains with A = 0, f(E) can be considered a 
constant, f(E) — /o, as long as \(E — n)/n\ <C 1. On 
the other hand, when A is finite, coherence peaks arise 
in the density of states^, so that at \E\ > A 



f(E) = 

fo 



E 



A 2 



(19) 



where in this expression E is measured relative to the 
center of the superconducting gap. [Eq. ( |19[ ) ignores the 
potential effect of thermal broadening of the coherence 
peaks.] 



The expression of Eq. (191 indicates that the conduc- 



-4 -2 2 4 

electron energy, e 



tance into a single grain is greatly enhanced when the 
voltage is aligned with the edge of the superconducting 
gap. For an array of grains, the total conductance is 
the integrated conductance of all the individual grains, 
each of which has a different relative alignment with the 
voltage. Thus, the differential conductance satisfies 

eV 

G(A,V) = G AD 1+d f 9l (E)f(eV-E + A)dE, (20) 



FIG. 6. (Color online) Single electron and pair DOGS, g*{e) 
and (/2(e), °f a 3d array of monodisperse metallic grains. The 
DOGS curves obey the same scaling relations as in 2d, [see Eq. 
( IT l] , indicating the presence of the same effective charges le, 
V2e, and 2e and A* = 0, A* = 1, and A* > 1, respectively. 



the DOGS can be observed directly from tunneling ex- 
periments. 

Tunneling experiments have previously been used 
to directly ob serve the Coulomb gap in lightly-doped 
semiconductors^ 2 ^ and have also measured the super- 
conducting gap in isolated superconducting grains^ and 
in disordered filma^. It is therefore natural to think that 
the single-electron DOGS gi(E) predicted here can also 
be measured via tunneling. In the problem we are con- 
sidering, however, the energy scales E c and A are similar 
in magnitude, and thus the tunneling conductance re- 
flects a convolution of the DOGS gi(E) with the density 
of states f(E) within each grain. As a result, we consider 
it worthwhile to explicitly state our predictions for the 
tunneling conductance G(A,V), where V is the applied 
voltage, at different values of the gap A. 

For simplicity, in this section we ignore the potential 
effects of spin polarization on the tunneling rates. This 
is equivalent to assuming that any applied magnetic field 
modifies the superconducting gap A primarily through 
orbital effects rather than the Zeeman effect, so that the 
electron energy levels shown in Fig. [2] are not labeled by 



where Go is a constant and A is the area of the tunnel 
barrier. [The term +A in the argument of / in Eq. 
accounts for the fact that the function f(E) in Eq. 



is defined relative to the center of the superconducting 
gap while the ground state energies described by gi{E) 
include the gap energy A.] 

Given our results for g%(E), one can use Eq. (20) to 
numerically evaluate the conductance G(A, V). For the 
limiting case A = 0, where the density of states f(E) 
is constant, Eq. (20) becomes simply G(A = 0, V) oc 

Iq V 9i(E)dE, or in other words g x {eV) cx dG(0,V)/dV. 
For small but finite A, on the other hand, such that < 
A* < 1, the conductance G(A, V) is enhanced at small 
V relative to G(0, V) as a result of the coherence peaks. 
At large A* > 1, a gap opens in gi(E), and G(A,V) 
remains at zero for \eV\ < (A* — 1)E C . 

This result is shown in Fig. [7] for the case of tunnel- 
ing into a 2d array. Here the conductance G(A, V) is 
plotted normalized to the value G(0, V) as a function of 
dimensionless voltage eV/E c for different values of A*. 
One can think that these different curves correspond to 
different magnetic field, since, as explained above, an in- 
creased magnetic field reduces the gap A. Thus, Fig. 
[7] suggests that if one starts with a sample for which 
A* > 1 and increases the magnetic field, a dramatic 
change occurs in the quantity G(A, V)/G(0, V). Namely, 
G(A, V)/G(0, V) first remains at zero for small V, since 
the single-electron DOGS is gapped. As the magnetic 
field is increased, the width of this gap decreases, un- 
til the point where A* = 1 and it disappears. Once 
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A* < 1, the value of G(A,V)/G(0,V) undergoes an 
abrupt change such that it becomes divergently large at 
small voltage. This divergence can be seen as the result 
of the coherence peaks, which greatly increase the tun- 
neling at small voltage relative to the case where there is 
no Cooper pairing (high magnetic field). Increasing the 
magnetic field also has the effect of lowering the conduc- 
tance peak at larger voltage, eV/E c ~ 2 + A*. 




2 3 4 

voltage, eV/E 



FIG. 7. (Color online) Tunneling conductance G(A, V) into 
a 2d array of superconducting grains as a function of volt- 
age V and for different values of the superconducting gap A*, 
plotted as a ratio of the conductance at A = 0. The volt- 
age V is plotted in the dimensionless form eV/E c . As the 
magnetic field is increased, driving down the value of A, the 
ratio G(A, V)/G(0, V) at small voltage changes from zero to 
a divergently large value as A* is made smaller than 1. 



It should be noted that the prediction of Fig. [jj is 
dependent on the existence of unscreened, long-range 
Coulomb interactions between electrons, which create the 
Coulomb gap in g\(E). If such long-range interactions 
are screened by the presence of a nearby metal electrode, 
which creates an image charge for each charged grain and 
truncates the 1 jr interaction, then the Coulomb gap will 
not be preserved and the predictions of this paper will 
be modified. For 2d arrays, it is therefore likely that 
macroscopic tunneling experiments will not be effective 
in identifying a Coulomb gap. The behavior of Fig. [7] 
may nevertheless still be identified if one uses a scanning 
tunneling tip to measure the conductance through indi- 
vidual grains and takes an ensemble average (as in, e.g., 
Ref . I5U|) . Alternatively, one can measure the conductance 
into a 2d face of a thick, 3d array, as was done in Ref. [2^1 
While Fig. [7] plots the conductance assuming tunneling 
into a 2d array, if one assumes that g\ (e) is identical to 
that of a 3d system (Fig. [6]) , the results are qualitatively 
very similar. 



VI. CONCLUSIONS 

In this paper we have proposed a model of a disordered 
granular superconductor and evaluated the DOGS and 
resistivity at different values of the superconducting gap 
A. Our primary result is the DOGS for single electrons 
and electron pairs shown in Figs, [l] and [6j We also have 
considered the implications of the DOGS for the conduc- 
tivity of the system (Figs. [3]), and explained a mechanism 
for negative magnetoresistance (Fig. |4| . Our predictions 
for the tunneling conductance are given in Fig. [JJ 

Perhaps the most remarkable result is the existence of 
effective charges le, 2e, and y/2e at A* = 0, A* > 1, and 
A* = 1, respectively, which was first reported by Ref. 
H5l These effective charges codify exact scaling relations 
between different results for the DOGS [Eq. (17)] and 
for the conductivity at low temperature [Eq. (18)], and 
can be understood in a fairly intuitive way. At A* = 0, 
electrons are unpaired and electronic conduction is per- 
formed by single electrons. At A* > 1, electrons become 
bound in Cooper pairs and these pairs are the primary 
players both in the conductivity and in determining the 
DOGS. At the point A* = 1, single electrons hop in a dis- 
order potential that is shaped primarily by Cooper pairs, 
and the single-electron DOGS and conductivity can be 
described by an effective charge e* = v2e. 

It is perhaps worth emphasizing that this effective 
charge y/2e does not represent a real quasiparticle in 
the traditional sense. For example, unlike the charges 
le and 2e, the charge v2e is unlikely to appear in the 
shot noise of the current (or the Fano factor) , since the 
actual hopping is performed by single electrons. Rather, 
the appearance of the charge \/2e in gi(E) and Tes is 
the result of a degeneracy in the electronic spectrum, 
which results in electrons being paired in the ground 
state. These paired electrons rearrange in the presence of 
a disorder potential and determine the properties of the 
ground state, while transport is carried by singles. It is 
this combination of intimidation by pairs and conduction 
by singles that produces the appearance of a \Jl charge. 

More generally, this view represents something of a 
novel paradigm in hopping transport. Namely, that a 
system can be simultaneously populated by two or more 
charged species (here, singles and pairs), one of which 
determines the Coulomb landscape while the other is re- 
sponsible for transport. Exploring this kind of physics in 
other classes of disordered systems remains a promising 
topic for future study. 
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